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Abstract 



■ We investigate e~^e ^ tt by taking into account wave function distortion ef- 

in 

Q^ ' fects in the near-threshold region and interpolating the correction factor to the 

Q^' high-energy perturbative QCD region. The strong color attraction between 



the color singlet t and t enhances the cross section over the tree level results 
near threshold. The rise of the cross section near the threshold is smoothed 



^ ■ by the effects of the large top width and below-threshold resonances. The 



reliability of the prediction is well assessed by comparing with experimental 
R ratios for lighter quark production at lower energies. The cross sections 
obtained here using the reported values of the top mass, can be used directly 
for New Linear Collider (NLC) projections. 
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I. INTRODUCTION 



Recent reports of the observation of the top quark, with a mass of 174 GeV [|T|, provoke 
improved theoretical studies of the top flavor. The planned and proposed experimental 
focus for the next decade to develop and understand top physics and top characteristics 
add impetus to theoretical investigations [0-0] • The super-heavy value of the top mass also 
opens up new vistas of electroweak physics and possibilities for probing for effects beyond the 
Standard Model. Many special features of the top have been studied recently PH^. Being 
heavier than the W, the top is extremely susceptible to weak decay by the t W^b channel. 
This large width distinguishes it from its doublet partner b, as well as the lower generation 
quarks and will lead to the disappearance of the sharp resonance structure corresponding 
to below-threshold resonances, as one transits from the c and b to the t flavor. Yet another 
feature, special to the top sector is the reliability of perturbative QCD to probe its bound 
and continuum states as the non-perturbative regime does not contribute effectively at the 
very short distances characterized by l/rrit for top flavor production. 

Careful analyses of higher-order QCD corrections for production of heavy quarks with 
a finite mass have been carried out [[TSHlSll which indicate that the next-to-leading order 
corrections are very large when the relative velocity between the top quark and the top 
antiquark is small. For such massive fermions near the threshold, the distortion of the 
outgoing relative wavefunctions between the top quark and the top antiquark, due to the 
strong color field connecting them, would be a dominant correction to the ti production 
processes. This can be understood in terms of 'Coulomb gluons' analogous to those acting 
between c and c in the J/tp systems and between b and 6 in T. The importance of such 
Coulomb effects for heavy quark production processes has been recognized earlier in the 
study of heavy quark productions in bound and continuum states |lT6|-[19|,p| . 

In the present work we expand on these aspects and predict the rates for ti production by 
virtual 7 and Z exchange in e~^e~ collider experiments, according to e^e~ 7*(Z*) QQ, 
including the important near-threshold region. The two loop modified running coupling 
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constant ots{<f) in the five flavor theory is used and the threshold QCD corrections are 
folded into the high energy results by an interpolation formula. Effects of top width and 
below-threshold resonances are included. Corrections for running of the QED coupling have 
been incorporated. The cross section obtained here using the latest top mass, can be used 
for future e^e~ collider projections. Results are extended to the charm and beauty flavors 
and applied to the familiar R parameter for comparison with experiment. 



II. TREE LEVEL CROSS SECTION 



The differential cross section for e+e QQ above threshold at a center-of-mass energy 



^ is given at the tree level by Eqs. (24. 11)- (24. 19) of Ref. [|0 



rfa(.) _ SvraVeQ y^^^ ^ (1) 



dz 4s V e 

where z = cos 6, and niQ and eg are the quark mass and charge, and we have included the 
three different colors of the quark- antiquark pair. In Eq. (|I]), the quantity j3 is 



/3=1/1 

the quantity D is 



4m 
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D = elA- 2eQXi[VVfA - 2af(3z] + Xi'B , (2) 



where 



A = l + z' + {l-l3'){l-z'), 
B = V^{1 + V^)A + (3'^a}{l + V^)[l + z"^] - 8(3VVfafZ, 

X 

16 sin^ 6yj cos^ 9^, ' 
V = -1 + 4sin2^^, 

l^/ = 2T3/-4eQsin2^^, 
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ttf = 2T3/ = 1 for the top, and 6^ is the Weinberg angle. The quantity x contains the Z 
mass mz as 



s 

X = 



s — m| 



Integrating over the angle, as we are not probing the asymmetry but the top production, 
we obtain 

(3) 

where 

H=(l + [el - 2eQXiVVj + xlVf{l + V^)]+ (5^a){\ + T/^) . (4) 

The above expressions are obtained using plane waves for all the fermions. It is however 
well recognized that the lowest order diagrams are inadequate to describe the quark sector, 
because of the large value of the QCD coupling constant as- Even for light quarks, the next- 
to-leading order radiative corrections for such production vertices predict an increase in the 
lowest order values, according to the correction factor (1 + cts/vr), where is the strong 



coupling constant |21|. These corrections, obtained painstakingly as the residual quantities 
after cancellation of singularities, unfortunately require the massless or chiral limit for their 
applicability. It is easy to see that the massless limit is not a desirable modus operandi when 
dealing the top flavor near the threshold, with its mass of 174 GeV. 

The large mass of the top quark assures the importance of the part of the vertex correction 
due to the exchange of longitudinal gluons between the top quark and the top antiquark at 
the 7*tt, or Z*ti vertex, near the ti threshold at around 350 GeV. This exchange can be 
represented by a potential that acts between the t and i. 



III. VERTEX CORRECTION FACTOR 



If we consider direct materialization of a quark-antiquark pair by decay of an intermediate 
virtual photon for instance, the process is characterized by a region of very small relative 

4 



Q-Q separations, at a Q-Q distance of ~ a/ y^. The effective part of the Q-Q interaction in 
this region of small relative separation associated with the production process is the inverse- 
r, Coulomb-like term. The linear part of the potential, that serves to effect the binding and 
impose confinement, controls the large r behavior p2| , |23[] and needs not be considered in 
the collapsed spatial zone comprising the annihilation or production vertex. This is further 
augmented by the short lifetime of the t and i that prohibits their moving away beyond the 
perturbative range. 

Upon expressing the one-gluon exchange potential between the quark and the antiquark 
in the 1/r Coulomb-like form, one can obtain exact solution for the wavefunction ||2^ and 
extract the r = part to correct the plane waves used in the Born diagram for their 
distortion in the color field. This procedure corrects for the distortion effect to all orders, 
and the resulting corrections can extend the applicability of PQCD for this process to the 
near-threshold region. We have discussed these correction factors for the lighter quark 



annihilation and production in a separate work ||25|. Though we restrict the main analysis 
to top flavor, the results can be easily extended to lighter flavors, and are seen below to be 
compatible with charm and beauty production near their relevant thresholds. 
We use a suitable short distance color potential of the Coulomb type, 

V{r) = 

r 

where 

is the appropriate color factor, and is the two-loop modified running QCD coupling 



constant and can be calculated from 



^ r,2^ ^ 47r / _ 2(51-19n^/3)ln[ln(gVA^)] \ 

^ (11 -2nj/3)ln(gVA2) V (11 - 2nj/3)2 In(g2/A2) )' ^ ' 

where is identified as s and A is taken to be A^^-* = 195 MeV for top production. 

For outgoing final states relevant to a produced QQ pair, the wave function of the quark 
Q in the field of the antiquark Q is given by pll 



tjj = iVe^P-'"(l - —ct ■ V)uF{-i^, 1, -i{pr +p-r)), 
2E 



(6) 



where u is the spinor for the free quark, F is the confluent hypergeometrical function and 
N is the normahzation constant given by 

27re 



1 - e-2-« 



where 



and V is the relative velocity between the quark and the antiquark given by ^5 



1 - iml/s 



(7) 



1 - 2ml/s ■ 

The square of the wavefunction taken at contact can be obtained as a generalization of 



the familiar 'Gamow-Sommerfeld factor' ||2^ as 

27r^ 



1 - e-2-€ 



(8) 



where spins of both quarks have been averaged over and the superscript (/) is to indicate 
that this is the correction factor for quark-antiquark production in the final state . 

In order to obtain a generalized correction factor, that is valid not only when Q and Q 
have low relative velocities, but also have large relative velocities, we use the interpolation 
technique suggested by Schwinger |28|. In the massless quark limit, the next-to-leading 
order QCD corrections to production of qq at high energies is given by the correction factor 
+ as/vr) E^. The transition of the QCD corrections from the correction factor 



of Eq. (j^) at low relative velocities, to the (1 + tts/vr) behavior at relativistic velocities 
can be accommodated with the introduction of an interpolative function f^-f\v), where 



cff 



- + v{-l + — 



(9) 



The generalized QCD correction factor, applicable at all relative velocities can be written 

as 



1 — exp[—2TTj(J)[v)\ 

Incorporating the distortion correction, equation @ for tt formation by the electroweak 
channels becomes, at energies above 2mQ, 



ais) = Kif^^-^m\l-'-^H (11) 



s \ e 



For virtual photon or Z mediated processes, the quark-antiquark pair must be in color- 
singlet states. The color factor Cf has the value 4/3 and the one-gluon exchange potential 
is attractive between the quark and antiquark. Due to this attractive nature of the color- 
singlet interaction, the correction factor (^) enhances top production over the tree level 
expectations. As in earlier analyses ||^, [|T^, there is a large enhancement of top flavor 
production near threshold due to the wavefunction distortion by this attractive color field. 



It may be noted that in the large ^ limit, the correction factor in Eq. (|TT|) takes the form 



TT, 



^ = nCjas/v. As i; — ^ 0, corrected a{s) of Eq. ([TT| ) does not diverge to infinity, because 



the Y 1 — 4mQ/s term of (0) cancels with the corresponding term in the velocity in Eq. (0). 
This cancellation of the Coulomb singularity by the phase space factor has been discussed 



in the context of QQ formation by qq annihilation in hadronic collisions [|T3[. In the absence 
of the vertex correction, the top production rate would plunge to zero at threshold, whereas 
with the correct inclusion of distortion effects, the cross section attains a relatively constant 
value. The introduction of the bound toponium states and the finite top width modifies 
this constancy as one approaches the threshold and softens considerably the cut-off at the 
threshold. 



IV. INCLUSION OF THE TOP DECAY WIDTH AND BELOW-THRESHOLD 

BOUND STATES 

The large weak decay rate of the top and the contribution of below-threshold bound 
states should be included for correct prediction of near-threshold top production. In the 
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Standard Model, the decay rate for the top by the dominant t W^h mode can be written 



as [E|] 



where 

L = [1 - {m^jmtff + [1 + ("^feMt)^][("^«>Mt)^] - 2{m^/mtY , 
( [mf - {rriy, + mb)'^Y^'^[mf - [m^ - mb) 

and V{th) is the Standard Model mixing parameter which can be taken to be unity. 

The effects of the top width can be introduced by including the contribution from all 
other energies to the cross section at the energy under consideration. The cross section at a 
given energy a/s should then be written as 

o{s) = Y.a,{Ei)Us,Ei), (13) 

i 

where ai{Ei) is the total cross section for state Ei and Ji{s,Ei) is the normalized Breit- 
Wigner distribution 

1 r/2 
^ ^ [(v/i - 2m, - E,y + (r/2)2] ■ 

The summation runs over all state energies E^ and the width F refers to the decay of the 
combined ti system and is given by F = 2Ft. The cross section can be devolved into a 
summation ab{s) over all discrete bound states and an integration o"c(s) over all continuum 
states as 

a{s) = (Tfe(s) + ads) . (14) 

Integrating over the continuum, the open or continuum cross section, including the top width 
effects can be written as 

,_^^s)=K'f>'-^(?fpHMs), (15) 
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Jc(s) = — i — h tan ' 



where 

.TJs) = , , , 

TT 12 V r 

In this case the width softens the sharp step-function type onset of continuum production 
into a smooth leakage into the region below the threshold. 

The case of the bound states requires the input of the bound state energies and wave- 
functions. It is expected that the attractive color field between the produced tt will form 
bound toponium states like the charmonium and bottonium systems with similar spectro- 
scopic properties. The tt system differs from the cc and bb in two important respects. Due to 
the large virtuality of the process for the pair with top flavor and the correspondingly large 
QCD scale, characterized by l/rrit, the QCD perturbative coupling is considerably smaller 
than that for the lighter qq pairs and PQCD is justified. Secondly, the large top mass and 
resulting large decay probability into the W and 6 is a characteristic not shared by the pairs 
with 6 or c flavor. 

The large magnitude of the momentum transfer allows us to describe the lowlying to- 
ponium bound states as those in the attractive perturbative Coulomb-like inter-quark color 
field as nonperturbative effects are not manifested at the short distances applicable to tt 
materialization. The effects of the gluon condensate have been studied by pr2|,|8| and shown 



to be small. The work in further establishes the validity of the Coulomb color potential 
for the toponia states. 

Bound states in the ti color field can be described simply by Coulomb-like wavefunctions 
and energy levels. Using the same perturbative Coulomb-like potential as for the continuum, 
we can write for the energy of the nth bound state 



2^2 ' 

where /i is the reduced mass of the quarks. The wavefunction at ti contact can be written 

m 



as 



TX \ n / 
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The total cross section for below-threshold resonance toponium production can then be 



written as 



" s \ e J An^ 
Introducing the resonances and the finite width, the cross section for the nth bound state 
becomes an{En)Jn{s, En) and F continues to be twice the top weak decay width as branching 
ratios to any other channel is negligible. The effect of the large magnitude of the top width 
of 1.5 GeV, is to smear out the resonance structure of the bound states and yield a smooth 
contribution. 

The total production cross section sums the contribution of bound and continuum states. 
The resulting continuum cross sections are enhanced over the tree level values, the effect 
being most pronounced just above threshold and tapering off at very high energies. This is 
demonstrated in Fig. 1, where the total cross section for e^e^ tt, calculated by including 
the width effects and the QCD correction factor A'*^-^^ for the continuum states, is shown 
as the solid curve and the tree-level cross section for the continuum states is shown as the 
dashed curve for comparison. The large enhancement at low relative tt velocities induced by 
the distortion corrections in this region is clearly visible. The cross sections remain enhanced 
over the lowest order results as the energy increases. At higher energies the cross sections 
fall with the usual inverse s behavior. 

We also present in Fig. 2 the detailed threshold behavior, with the separate contributions 
from the bound and the continuum states. As remarked earlier, the large top width smears 
out the resonance structure into a smooth contribution. These results are in conformity to 
earlier investigations of the threshold region, that were spread over wide-ranging possible 
top mass values [3^TT||. 

We have used the reported central value of 174 GeV for the top mass [|l|. This value is 
in line with the earlier direct searches that set a lower limit > 131 GeV/c^ and with 
the global fits to precision electroweak fits that give a central value of rrit = 177 GeV/c^ |^ 



The accepted value of 91.17 GeV has been used for the Z mass and the weak mixing 
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angle has been taken as rj = 0.225 p2|,p6 . 



The correction factor K^'^^ used here is similar to those used earlier p!6|-p!9|,p|,p8| , for lighter 
quarks, with small differences. We use relativistic kinematics to interpolate the correction 
factor to the high-energy region. Our results are in general agreement with the earlier works, 
but serve to quantify and extend them using the latest values of the electroweak parameters 
and top mass and are directly applicable over a wide range of energies. 

V. LIGHTER HEAVY QUARKS 

In order to demonstrate the reliability of our predictions for top quark production, it 
is necessary to show that the large enhancement of heavy quark production is consistent 
with experimental data. Therefore, we extend our analysis to cover the lower energy region 
spanned by experiments for heavy quark formation by e~^e~ annihilation. 

The most appropriate quantity to use for meaningful comparison with experiment is the 
ratio of hadronic to leptonic (muonic) branching ratios, commonly designated as R, 

^ cr(e+e^ — > hadrons) 
cr(e+e~ — > f^~^f^~) 

The numerator takes the form 

a{e+e- hadrons) = N,J2 K'-f\qf)af{s), 

f 

where for the continuum states 



In this case the total widths being of Mev order have no perceptable influence on flavor 
production cross sections. Nor do they distort the resonance levels, which can be introduced 
at the appropriate below-threshold points. 

The results of the R ratio calculated for the continuum states by including the correction 
factor K^-l"^ using A = A*^^) = 338 MeV in Eq. for c and b production are shown as 
the solid curve in Fig. 3 and are compared with the experimental data The R ratio 
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obtained just from the tree-level diagram without the correction factor is also shown as the 
dashed curve. The results of Fig. 3 show that the inclusion of the correction factor gives a 
good description of the general features of the R ratio. The experimental R ratio and the 
corrected ratio R are much greater than the tree-level results, near the threshold region. The 
large enhancement of heavy-quark production is consistent with experimental data, lending 
support to the predictions of Fig. 1 and Fig. 2 for ti production. 

The exhaustive fits to the energy dependence of R beyond charm threshold by Barnett 
et al [|l7l also used the interpolation formula due to Schwinger and included the low energy 
dependence of the cross section through an inverse velocity term. They used smearing 
procedures for proper match between theory and experimental data above the threshold and 
investigated new hypothetical particles. They expect their threshold results to be unreliable 
due to their inclusion of the distortion effect only up to the first order in a^. In the present 
work, we incorporate the distortion corrections to all orders of the coupling constant on 
the one hand and we interpolate to the massless quark limit on the other hand. The 
simple analytical semi-empirical correction factor K^-^^ of Eq. (^) obtained with only a 
Coulomb-like potential reproduces the gross features of R over a coarse mesh of energy, 
as the energy scale of Fig. 3 and the need for the smearing procedure of Ref. will 



indicate. In finer energy resolutions, the experimental R ratio is modulated by the presence 
of resonances above the continuum threshold which necessitates the description with a more 
refined potential including the additional influence of the long-range part of Q-Q interaction, 
as discussed well in Refs. p2| , |33| , |3^ . The agreement of our results with the gross features 



of the experimental R ratio adds credence to the approximate validity of the calculated tt 
production results. 

VI. DISCUSSION. 

The importance of utilizing the largest possible cross section for top production is self 
evident. In the hght of possible manifestations of new physics expected for the top era 
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and Standard Model violating effects [^], it is important to have a meaningful sample of top 
events. It is also imperative to have a clear quantifications of Standard Model predictions 
to enable clean extraction of exotic physics signals. Our distortion corrected cross sections 
provide the expected top yield within the framework of the Standard Model and can serve 
as the correct reference to probe new physics frontiers that might be discernible beyond the 
top threshold. 

The large top mass introduces another feature of interest that merits discussion. This is 
the coupling of the Higgs boson to the ti pair, that is manifested as exchange of the scalar 
Higgs boson between the members of the produced pair. Such an exchange leads to an 
additional attractive potential between the t and t that takes the Yukawa form and can be 



written 35 



1 /mA^e""""'' 

VH{r) 



47r V V / r 

where mn and v are the mass and the vacuum expectation value of the Higgs. This is 
expected to enhance the attractive effects of gluon exchange in the color singlet combination. 
However in the absence of reliable predictions for the Higgs mass, inclusion of these effects 
may be premature. Projections for Higgs effects have been carried out by 0. 

The possibility of toponium states bound in the pure Higgs potential could also arise. 
Bounds for the exchanged particle mass, to allow for bound states in a Yukawa potential 
can be estimated [^^. The present lower limit of 60 GeV for the Higgs mass |^ rule out 
any exclusive bound states in the Higgs potential for t and t of the detected mass range. 

Finally a comment on possible variations in the top width seems called for. It is well 
known that initial state binding and masses of the decay products are connected to the 



effective decay rate of a fermion ||37| , |38| . Since in this case, the binding energy of the 
decaying top to its production partner, i, is of the order of 0.1 GeV only, these effects 
are small and expected to be beyond the sensitivity range of experimental access at this 
time. The influence of momentum dependence on the top width and production has been 
addressed recently p9|. 



13 



It may be pointed out that the enhancement due to the exclusive color-singlet tt formation 
in electroweak production will not occur in qq annihilation processes in pp and pp collisions 
at Tevatron and LHC Nevertheless, the influence of distortion effects in these environments 
merits a discussion. The gluon mediated qq QQ process is suppressed by the distortion 
correction E^ due to the repulsive color-octet nature of the one-gluon exchange potential 
between q and q at the initial vertex and between Q and Q in the final vertex. However the 
magnitude of this suppression is much smaller than the enhancement for the color-singlet 
electroweak channels because of the different color factors for the two cases. The distortion 
effects are not as striking for QCD vector boson mediated qq processes as for electroweak 
boson mediated ones. However in hadronic and nuclear collisions, qq processes are in any case 
suppressed compared to gluon induced ones. We are investigating the distortion corrections 
on the gluon fusion channels in a separate work. 

We reiterate the justification of the use of the perturbative color potential as acting 
between both bound and continuum states due to the double influence of the large virtuality 
of the process, corresponding to contributing distances of order l/m^ and the large top decay 
rate which serves to kill the top before it can dress itself or cover the wider range of toponia 
levels. Our results, using the experimentally reported top mass can be used directly for 
projections of top flavored event rates at the proposed New Linear Collider (NLC) facility 
and may help to select the most efficient colliding energies for productive top hunting by 
e~^e~ annihilation. 
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FIGURES 

FIG. 1. The cross section for e^e~ — > ti. The sohd curve is the total cross section by including 
the bound and continuum states, and the correction factor K^^\ The dashed curve is obtained 
from the tree-level diagram for the continuum states. 

FIG. 2. The cross section for e"'"e~ near the threshold region. The solid curve is the total 
cross section by including the bound states and the continuum states. The dot-dashed curve and 
the long-dashed curve are obtained from the bound states and the continuum states respectively. 

FIG. 3. The ratio i? as a function of the center-of-mass energy ^/s. The solid curve is the 
result obtained by including the correction factor K^^'> , and the dashed curve is from the tree-level 



diagram. The data points are from the compilation of Ref. |26 
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